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Abstract. With a dynamical mass Md yn ~ 1.3 x 10 5 Mq and a lower limit M c \ > 5 x 10 4 Mq 
from star counts, Westerlund 1 is the most massive young open cluster known in the Galaxy and 
thus the perfect laboratory to study massive star evolution. We have developed a comprehensive 
spectral classification scheme for supergiants based on features in the 6000 - 9000A range, which 
allows us to identify > 30 very luminous supergiants in Westerlund 1 and ~ 100 other less evolved 
massive stars, which join the large population of Wolf-Rayet stars already known. Though 
detailed studies of these stars are still pending, preliminary rough estimates suggest that the 
stars we see are evolving to the red part of the HR diagram at approximately constant luminosity. 
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1. Introduction 

As they evolve towards the Wolf-Rayet (WR) phase, massive stars must shed most 
of their outer layers. Models predict that massive stars will become supergiants (SGs) 
and evolve redwards at approximately constant Lboi, but observations reveal a complex 
zoo of transitional objects, comprising Blue SGs, Red SGs, Yellow Hypergiants (YHGs), 
Luminous Blue Variables (LBVs) and OBfpe/WNVL stars, whose identification with 
any particular evolutionary phase is difficult. Understanding this evolution is, however, 
crucial because the mass loss during this phase completely determines the contribution 
that the star will make to the chemistry of the ISM and even the sort of post-supernova 
remnant it will leave. 

Unfortunately, massive stars are scarce and, as this phase is very short on evolutionary 
terms, examples of massive stars in transition are rare. For most of them, distances are 
unknown and so luminosities arc known at best to order-of-magnitude accuracy, resulting 
in huge uncertainties in other parameters (M», R+). The difficulty to place these objects 
in the evolutionary sequence is obvious. 

The open cluster Westerlund 1 (Wd 1) offers an unprecedented opportunity to improve 
this situation. It contains a large population of evolved stars which have formed at the 
same time, are at the same distance and have the same chemical composition. With an 
age of 4-5 Myr, Wd 1 contains a rich population of stars in transitional states that can 
be used to constrain evolutionary models. 

f Present address: Chester F. Carlson Center for Imaging Science, Rochester Institute of 
Technology, 54 Lomb Memorial Drive, Rochester NY, 14623, USA 
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Figure 1. Sequence of /-band spectra of blue supergiants in Wd 1, showing the evolution of the 
main features. Note the disappearance of N I features between Pa 12 and Pa 13 and O I 8446A 
around B2 and the development of Pa 16, which shows the presence of a strong C ill line for 
stars B0 and earlier.. 



2. Cluster parameters 

The parameters of Wd 1 are still poorly determined, though significant progress 
has been made in recent years. Two estimates of its mass have been recently made. 



Mengel & Tacconi-Garman (2007) have used the radial velocity dispersion of the ten 
stars brightest in the infrared (a = 8.4 km s -1 ) to estimate a mass of ~ 1.3 x 10 5 Mq. 
This estimate suffers from a large uncertainty due to the low number of stars used and 
the spectral peculiarities of some. Brandner et al. (2008) have used star counts in the 
infrared to set a lower limit on the cluster mass M c \ > 5 x 10 Mq. Again, there are im- 
portant uncertainties coming into this estimation. Large discrepancies have been found in 
the determinations of distance and (hence) masses based on infrared pre-main-sequence 
tracks and optical post-main-sequence tracks. In any case, the integrated initial mass 
function appears consistent with Salpeter's and then the observed population of massive 
stars would imply a mass ~ 10 5 Mq, consistent with both determination. 

An important source of uncertainty in both measurements is the role of binarity, which 
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affects both radial velocities and transformation of magnitudes into masses. Recent re- 
sults suggest that the fraction of binaries with similar mass components is very high 
amongst cluster members. Most Wolf-Rayet stars show indirect indications of binarity, 
presence of dust in WC stars (Crowther et al. 2006) and hard X-ray emission in WN 
stars (Clark et al. 2008) Some of them also show photometric variability and at least 



one is an eclipsing binary (Bonanos 2007) The identification of a large number of X-ray 



sources detected by Chandra with a population of evolved late-0 stars (Clark et al. 2008) 
suggests that the high binary fraction extends to lower masses. 

The distance to the cluster is not very well constrained either. Photometry is affected 
by the very strong reddening (Ay « 12 mag). Analysis of the E(B — V) colours suggests 
that the reddening deviates from the standard law (Clark et al. 2005) A recent determi- 
nation, making use of atomic hydrogen in the direction to the cluster, gives d — 3.9 ± 0.7 
(Kothes fc Dougherty 2007)| compatible with, though slightly shorter than, estimates 
based on the stellar population (e.g., Crowther et al. 2006). 

The age of the cluster can be constrained from the observed population. The ra- 
tio of WR stars to red and yellow hypergiants favours an age of 4.5 or 5.0 Myr, with 
the progenitors of the Wolf-Rayet stars having initial masses in the 40 - 55 Mq range 
(Crowther et al. 2006) Such age is fully compatible with the observed population of 



blue supergiants (see below), which should be descended from stars with initial masses 
of ~ 35 Mq. This again would imply masses of ^ 30 Mq for the stars at the top of the 
main sequence, in good agreement with the 07-8 V spectral type, again appropriate for 
the age. 



3. Observations and analysis 

Observations of stars in Wd 1 were carried out on the nights of 2003 June 12th and 13th 
using the spectro-imager FORS2 on Unit 1 of the VLT (Antu) using three different modes: 
longslit, multi-object spectroscopy with masks (MXU) and multi-object spectroscopy 
with movable slitlets (MOS). We used grisms G1200R and G1028z to obtain intermediate 
resolution spectroscopy. With this setup, we obtained almost continuous coverage over 
the 5800-9500A range at intermediate resolution. 

Within the central 5' x 5' field of view, we selected our targets from the list of likely 
members of Clark et al. (2005). For the external regions, targets were selected at random 
amongst relatively bright stars. In total, we took three MXU and one MOS mask with 
both G1200R and G1028z grisms, and two further MXU masks with only the G1200R 
(these were aimed at relatively faint objects, which were expected to be OB stars near 
the MS and so not to have strong features in the range covered by G1028z). This resulted 
in ~ 100 stars observed with G1200R and ~ 70 stars observed with G1028z. More than 
90% of the spectra turned out to correspond to OB stars, and hence cluster members. 

We have used these observations to derive spectral types for cluster members. The use 
of /-band spectra to classify OB stars has been explored by Caron et al. (2003) and is 
further discussed in Appendix A of Clark et al. (2005). We have studied further spectral 
type and luminosity indicators in the R and / bands and derived spectral types from 
a combination of features. Features that contribute to our analysis are: the shape and 
strength of emission (P-Cygni profiles) in Ha, the strength of the Cn 6578,6582A doublet 
(in absorption), the presence and strength of the Cn 7231, 7236 A doublet (in emission, 
a wind feature), the strength of the O I 7774 A triplet and the strength of the Nil 6482 A 
line. In the / band, apart from the shape and strength of the Paschen lines (see Clark et 
al. 2005), a main indicator is the presence of the C in 8502 A line, which appears, blended 
into Pa 16, for stars B0 and earlier (see Fig. [I]). 
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Figure 2. Sequence of red spectra of blue supergiants in Wd 1. The shape of Ha changes 
from a P-Cygni profile to pure emission around spectral type BO. 5. The N II 6482A line is only 
prominent in la supergiants and its strength peaks sharply at B2. 



4. Implications 

The upper HR diagram. The stellar content revealed by these observations is as fol- 
lows. There is an A supergiant showing spectral variability (W243; Clark & Negueruela 
2004) and six other A/F hypergiants of very high luminosity, with My « —9 to —10 



(Clark et al. 2005) We call these objects hypergiants because of their very high intrinsic 



magnitudes, without making any direct inference of their evolutionary status. 

The same can be said of at least three late-B luminous stars: W42a (B8Ia + ), W33 
(B5Ia + ) and W7 (B5Ia + ). In addition, there are ~ 20 supergiants that can be unam- 
biguously classified as la, covering the B0Ia-B4Ia range. The earliest very luminous 
star that we find is W74, with spectral type 09.5 and a luminosity class that could be 
la or lab. 

Even though the exact luminosity class of O-type supergiants is more difficult to deter- 
mine using red spectra, we do not find evidence for any other luminous O-type supergiant, 
though there are > 50 and (taking into account the incompleteness of our dataset) likely 
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~ 100 stars which we generically classify as 091, meaning late-0 stars well above the 
main sequence. 
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Figure 3. Spectra of three transitional objects in Westerlund 1, which can be classified as very 
late WN stars or early-B hypergiants. These objects might be examples of stars looping back to 
the blue region from the red supergiant phase. 



Evolution . We have made a rough attempt at calculating the intrinsic luminosities 
of all the supergiants, by using their V and / magnitudes and assuming that the re- 
lation between E(V — I) and Mi is standard. We use the intrinsic (V — I) colours of 



Ducati et al. (2001) and bolometric corrections from Martins et al. (2005) and Humphreys & McElroy (1984) 
Even though there is a large scatter, we find no clear trend. Stars of all spectral types 
seem to have similar bolometric magnitudes, a result consistent with the idea of red- 
wards evolution at constant bolometric luminosity. Note that we have not compared the 
bolometric magnitudes of the red supergiants, as their spectral types are uncertain, and 
the bolometric corrections of such bright red stars are unknown. 

This result does not imply that there are no stars looping back towards the blue in 
Wd 1. Obviously, the WR stars must have come back to the blue before reaching their 
present stage. There are three emission-line stars in the cluster which can be classified as 
very late WN stars or extreme B supergiants (see Fig- EJ) • These objects might be evolving 
into Wolf-Rayet stars. Detailed analysis of their spectra, in order to derive chemical 
abundances, will be necessary before their evolutionary status may be ascertained. Similar 
analyses should be conducted for the A and F hypergiants and W243. 

The next steps in our work involve an investigation of the binary fraction in the cluster 
and detailed analysis of the spectra with state-of-the art models in order to obtain a better 
constraint on their evolutionary status. 
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Discussion 

Lang: What does the distribution of radio emission reveal in this cluster - both the 
extended emission and the point-like sources? 

Negueruela: There are point sources associated with most of the cool and cold su- 
pergiants, which, we believe, are due to photoionisation of their extended atmospheres 
or winds by the hot stars. There is no extended emission obviously associated with the 
cluster, suggesting that supernova explosions have blown away any diffuse material. A 
more detailed summary may be found in the contribution by Sean Dougherty and Simon 
Clark to the proceedings of Massive Stars: Fundamental Parameters and Circumstellar 
Interactions, arXiv : 0705 . 0971vl. 

Crowther: As your title suggests, the presence of large numbers of massive stars in 
Wd 1 allows robust tests of evolutionary models. By way of example, N(WN with H) < 
N(WC) < N(WN without H), yet both single star (i.e., Geneva) and binary models (e.g., 
Eldridge, this conference) predict the completely wrong subtype distribution of N(WN 
without H) < N(WN with H) < N(WC) for an instantaneous burst of 4-5 Myr. 

Negueruela: Indeed. The observed population is most likely compatible with a single 
burst of star formation and hence offers an observational test of evolutionary models. 
Also, as I mentioned, the fact that we can observe a large population of massive (M» > 
30A/q) stars evolving away from the main sequence and a very large population of low 
mass stars evolving towards the main sequence gives us an excellent opportunity to try 
to set post and pre-main-sequence isochrones in the same reference time and compare 
their predictions. 



